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Abstract 



We have studied two-body charmless hadronic decays of B mesons into 
the final states mr, Kn, and KK. Using 3.3 million BB pairs collected 
with the CLEO-II detector, we have made the first observation of the 
decays B° — ► K + n~, B + — > K°tt + , and the sum of B + — > 7r + 7r° and 
B + — > K + tt° decays (an average over charge-conjugate states is always 
implied). We place upper limits on branching fractions for the remaining 
decay modes. 

PACS numbers:13.25.Hw,14.40.Nd 
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The phenomenon of CP violation, so far observed only in the neutral kaon system, can be 
accommodated by a complex phase in the Cabibbo-Kobayashi-Maskawa (CKM) quark-mixing ma- 
trix jij. Whether this phase is the correct, or only, source of CP violation awaits experimental 
confirmation. B meson decays, in particular charmless B meson decays, will play an important 
role in verifying this picture. 

The decay B° — > tt + it~, dominated by the b — ► u tree diagram (Fig. |](a)), can be used to 
measure CP violation due to B° - B° mixing at both asymmetric B factories and hadron colliders. 
However, theoretical uncertainties due to the presence of the b — > dg penguin diagram (Fig. 0(b)) 
make it difficult to extract the angle a of the unitarity triangle from B — > 7r + 7r~ alone. Additional 
measurements of B + — > 7r + 7r°, B° — ► 7r°7r°, and the use of isospin symmetry may resolve these 
uncertainties ||. 

B — > Kir decays are dominated by the b — > sg gluonic penguin diagram, with additional 
contributions from b — > u tree and color-allowed electroweak penguin (Fig. 0(d)) processes. In- 
terference between the penguin and spectator amplitudes can lead to direct CP violation, which 
would manifest itself as a rate asymmetry for decays of B and B mesons. Recently, the ratio 
R = B(B — > K ± tt t )/B(B ± —y K ^), was shown |J to constrain 7, the phase of V u b- Several 
methods of measuring 7 using only decay rates of B — > Ktt, tttt processes were also proposed [||. 
This is particularly important, as 7 is the least known parameter of the unitarity triangle and 
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FIG. 1. The dominant decay processes are expected to be (a) external W-emission, (b) gluonic 
penguin, (c) internal W-emission, (d) external electroweak penguin. 

is likely to remain the most difficult to determine experimentally. This Letter describes the first 
measurement of exclusive charmless hadronic B decays. Previous measurements existed only for 
the sum of several two-body final states [|]|| . 

The data set used in this analysis was collected with the CLEO-II detector at the Cornell 
Electron Storage Ring (CESR). It consists of 3.14 fb" 1 taken at the T(4S) (on -resonance) and 
1.62 fb" 1 taken below BB threshold. The on-resonance sample contains 3.3 million BB pairs. The 
below-threshold sample is used for continuum background studies. 

Charged tracks are required to pass track quality cuts based on the average hit residual and 
the impact parameters in both the r — (p and r — z planes. Pairs of tracks with vertices displaced 
by at least 3 mm from the primary interaction point are taken as K$ candidates. We require the 
7r + 7r" invariant mass to be within 10 MeV, two standard deviations (a), of the Kg mass. Isolated 
showers with energies greater than 30 MeV in the central region of the Csl calorimeter and greater 
than 50 MeV elsewhere, are defined to be photons. Pairs of photons with an invariant mass within 
20 MeV (~ 2a) of the nominal ir° mass are kinematically fitted with the mass constrained to the 
7T° mass. To reduce combinatoric backgrounds we require the lateral shapes of the showers to be 
consistent with those from photons, and that |cos#*| < 0.97, where 9* is the angle between the 
direction of flight of the ir° and the photons in the tt° rest frame. 

Charged particles are identified as kaons or pions using dE/dx. Electrons are rejected based on 
dE/dx and the ratio of the track momentum to the associated shower energy in the Csl calorimeter. 
We reject muons by requiring that the tracks do not penetrate the steel absorber to a depth greater 
than five nuclear interaction lengths. We have studied the dE/dx separation between kaons and 
pions for momenta p ~ 2.6 GeV/c in data using D* + -tagged D° — > K~ir + decays; we find a 
separation of (1.7 ± 0.1)cr. 

We calculate a beam-constrained B mass M = J ' — p 2 B , where pb is the B candidate mo- 
mentum and E^ is the beam energy. The resolution in M ranges from 2.5 to 3.0 MeV/c 2 , where 
the larger resolution corresponds to decay modes with 7r°'s. We define AE = E\ + E^ — E^, where 
Ei and E% are the energies of the daughters of the B meson candidate. The resolution on AE is 
mode-dependent and ranges from ±26 MeV for Kgir + to +82/ — 162 MeV for 7r°7r°. The latter 
resolution is asymmetric because of energy loss out of the back of the Csl crystals. The energy 
constraint also helps to distinguish between modes of the same topology. For example, AE for 
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TABLE I. Experimental results and theoretical predictions [10]. Branching fractions (23) and 



90% C.L. upper limits are given in 10 5 units. Quoted significance of the fit results is statistical 
only. The errors on B are statistical, fit systematics, and efficiency systematics respectively. 



Mode 


N s 


Sig. 


£{%) 


B 


Theory B 


7T + 7T~ 


9.9±|;S 


2.2a 


44 ±3 


< 1.5 


0.8-2.6 


vr+vr 


11.3±f;| 


2.8a 


37 ±3 


< 2.0 


0.4-2.0 


ttV 


2.71?-? 


2Aa 


29 ±3 


< 0.93 


0.006-0.1 


K+tt- 


21.6±g5 


5.6ct 


44 ±3 


1.5±g;l ± 0.1 ± 0.1 


0.7-2.4 


K+TT° 


8.711J 


2.7a 


37 ±3 


< 1.6 


0.3-1.3 


K°TT+ 


9.2±|;I 


3.2a 


12 ±1 


2.3±{- 1 ±0.3 ±0.2 


0.8-1.5 


K°TT° 




2.2a 


8± 1 


< 4.1 


0.3-0.8 


K + K- 


o o +1 - 3 

u - u -o.o 


0.0a 


44 ±3 


< 0.43 




K+K° 


u -°-0.6 


0.2a 


12 ±1 


< 2.1 


0.07-0.13 


K°K° 







5± 1 


< 1.7 


0.07-0.12 


h+TT° 


2o.oi|:| 


5.5a 


37 ±3 


1.61a! ± 0.3 ± 0.2 





B° — > K + ir , calculated assuming B° — > 7r + 7r , has a distribution that is centered at —42 MeV, 
giving a separation of 1.6c" between B° — > X + 7r _ and B° — > 7r + 7r~. We accept events with M 
within 5.2-5.3 GeV/c 2 and |A£?| < 200(300) MeV for decay modes without (with) a vr° in the final 
state. This fiducial region includes the signal region, and a sideband for background determination. 

We have studied backgrounds from b —* c decays and other 6 — ► u and 6 — > s decays and 
find that all are negligible for the analyses presented here. The main background arises from 
e + e~ — > (where q = u,d, s, c). Such events typically exhibit a two-jet structure and can produce 
high momentum back-to-back tracks in the fiducial region. To reduce contamination from these 
events, we calculate the angle 6s between the sphericity axis of the candidate tracks and showers 
and the sphericity axis of the rest of the event. The distribution of cos 9s is strongly peaked at ±1 
for qq events and is nearly flat for BB events. We require | costs') < 0.8 which eliminates 83% of 
the background. Using a detailed GEANT-based Monte-Carlo simulation g] we determine overall 
detection efficiencies (£) of 5 — 44%, as listed in Table ffl. Efficiencies contain branching fractions 
for K° — > K s — * 7r + 7r~ and 7r° — ► 77 where applicable. We estimate a systematic error on the 
efficiency using independent data samples. 

Additional discrimination between signal and qq background is provided by a Fisher discrimi- 
nant technique as described in detail in Ref. |J. The Fisher discriminant is a linear combination 
T = YU=i a i1li wnere the coefficients on are chosen to maximize the separation between the signal 
and background Monte-Carlo samples. The 11 inputs, yi, are | cos 9 can d\ (the cosine of the angle be- 
tween the candidate sphericity axis and beam axis), the ratio of Fox- Wolfram moments H2/H0 ||, 
and nine variables that measure the scalar sum of the momenta of tracks and showers from the 
rest of the event in nine angular bins, each of 10°, centered about the candidate's sphericity axis. 

For all modes except B° — ► K°K° we perform unbinned maximum-likelihood (ML) fits using 
AE, M, J 7 , I cos 6b\ (the angle between the B meson momentum and beam axis), and dE/dx 



5 



(where applicable) as input information for each candidate event to determine the signal yields. 
Five different fits are performed, one for each topology (h + h~ , h + ir°, 7r°7r°, h + Kg, and Kgir , 
referring to a charged kaon or pion) . In each of these fits the likelihood of the event is parameterized 
by the sum of probabilities for all relevant signal and background hypotheses, with relative weights 
determined by maximizing likelihood function (£). The probability of a particular hypothesis is 
calculated as a product of the probability density functions (PDF's) for each of the input variables. 
The PDF's of the input variables are parameterized by a Gaussian, a bifurcated Gaussian, or a sum 
of two bifurcated Gaussians, except for | cos#b| (1 — | cos#b| 2 for signal, constant for background), 
background AE (straight line), and background M (f(M) oc My/1 — x 2 exp[— 7(1 — x 2 )]; x = 
M/E h ) 0. 

The parameters for the PDF's are determined from independent data and high-statistics Monte- 
Carlo samples. We estimate a systematic error on the fitted yield by varying the PDF's used in 
the fit. The error is dominated by the limited statistics in the independent data samples we used 
to determine the PDF's. Further details about the likelihood fit can be found in Ref. ||. 

Figure || shows contour plots of — 21n£ for the ML fits to the signal yields (N). The curves 
represent the na contours (n = 1 — 5), which correspond to the increase in — 21n£ by n 2 . The 
dashed curve marks the 3<r contour. The statistical significance of a given signal yield is determined 
by repeating the fit with the signal yield fixed to be zero and recording the change in — 21n£. To 
further illustrate the fits, Fig. ^ shows M (AE) projections for events in a signal region defined by 
\AE\ < 2a ( \M — 5. 28 1 < 2aj\j). We also make a cut on T which keeps 67% of the signal and 
rejects 80% of the background. For Fig. ||(a), events are sorted by dE/dx according to the most 
likely hypothesis. For Fig. ^(c), 3<r consistency with the pion hypothesis is required. Overlaid on 
these plots are the projections of the PDF's used in the fit, normalized according to the fit results 
multiplied by the efficiency of the additional cuts (~ 60 — 70% for the signal and ~ 2 — 10% for the 
background). The central values of the signal yields from the fits (N$) are given in Table |. We 
find statistically significant signals for the decays B° — > K + ir~ and B + — * K®tt + . The latter mode 
constitutes the first unambiguous observation of a gluonic penguin decay. The former mode may 
have a sizeable contribution from the color-allowed b — > u tree-level spectator diagram in addition 
to the dominant gluonic penguin amplitude. We also observe a significant signal in the sum of 
decays B + K + ir° and B + -> vr+vr . 

As a cross-check, we perform a counting analysis in the modes B° — > K + n~, B + — > K°ir + , 
and B + — ► h + ir°. We calculate the probability of the background fluctuation to produce the excess 
of events shown in Fig. |3| to be 2.0 x 10 -7 for the K + tt~ mode, 1.6 x 10~ 3 for the h + ir mode, and 
2.5 x 10~ 4 for the K°ir+ mode. 

The statistical significance of the fitted yields in the modes it + it~, 7r + 7r°, 7r°7r°, K + ir°, and 
K°ir° ranges from 2.2a to 2.8<r. We consider these to be not statistically significant and calculate 
90% confidence level (C.L.) upper limit yields by integrating the likelihood function 

f NUL C^WdN 

IZ°£^(N)dN U ' yU [L) 

where jC max (N) is the maximum C at fixed N to conservatively account for possible correlations 
among the free parameters in the fit. We then increase upper limit yields by their systematic errors 
and reduce detection efficiencies by their systematic errors to calculate branching fraction upper 
limits given in Table I. 

We search for the decay B° — ► K°K° via K°,K° — ► Kg — ► ir + ir~. Since the background 
for this decay is quite low, the complication of a ML fit is not necessary and a simple counting 
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FIG. 2. Contours of the -21n£ for the ML fits to (a) N K ±^ and N n+n - for B° -> K + n~ 
and B° -» Tr+Ti-; (b) 2V K7r o and iV^o for 5+ -> and -► ^+7r°; (c) N K o K and JV^ for 

B + -> and 5+ ET°7r + . 
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5.22 5.24 5.26 5.28 -0.1 0.1 0.2 

M (GeV) A E (GeV) 

FIG. 3. M and AE plots for (a) B° -» K+tt-, (b) B + -» /i+vr , and (c) B + -► K°7r+. The 
scaled projection of the total likelihood fit (solid curve) and the continuum background component 
(dotted curve) are overlaid. 

analysis is used. Event selection is as described above, except no Fisher discriminant is used 
and |cos#t| < 0.75 cut is applied (cos#t is defined similar to cos#s, but with thrust axis used 
instead of sphericity). We define the signal region by requiring \AE\ < 65 MeV (2.5a), and 
\M — 5.28| < 0.005 GeV/c 2 (2.4<r). We observe no events in the signal region and calculate a 90% 
C.L. branching fraction upper limit of B(B° — ► K°K°) < 1.7 x 10~ 5 . 

As a comparison, we relate B — > nlv and B — > irir processes within the factorization hypothesis. 
Using the ISGW II || form factors, the QCD factor ax = 1.03 ± 0.07 ||, and the CLEO 
measurement B(B° -» k~1 + v) = (1.8 ± 0.4 ± 0.3 ± 0.2) x 1(T 4 @, we predict B(B° -> vr+vr-) = 
(1.2 ± 0.4) x 10" 5 and B{B + vr+vr ) = (0.6 ± 0.2) x 10" 5 These predictions are consistent 
with our upper limits as well as central values from the fit: B(B° — ► ir + ir~) = (0.7 ± 0.4) x 10 -5 
and B{B + -> tt+tt ) = (0.9±g;f) x 10~ 5 . 

In summary, we have measured branching fractions for two of the four exclusive B — ► Kir 
decays, while only upper limits could be established for the processes B — ► -kit, KK. Our results 
therefore indicate that the b — > sg penguin amplitude dominates charmless hadronic B decays. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent luminosity 
and running conditions. This work was supported by the National Science Foundation, the U.S. 
Department of Energy, the Heisenberg Foundation, the Alexander von Humboldt Stiftung, Research 
Corporation, the Natural Sciences and Engineering Research Council of Canada, and the A. P. Sloan 
Foundation. 
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